Low-lying dipole resonance in neutron-rich Ne isotopes 
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Microscopic structure of the low-lying isovector dipole excitation mode in neutron-rich 26 > 28 ' 30 Ne 
is investigated by performing deformed quasiparticle-random-phase-approximation (QRPA) calcula- 
tions. The particle-hole residual interaction is derived from a Skyrme force through a Landau-Migdal 
approximation. We have obtained the low-lying resonance in 26 Ne at around 8.5 MeV. It is found 
that the isovector dipole strength at E x < 10 MeV exhausts about 6.0% of the classical Thomas- 
Reiche-Kuhn dipole sum rule. This excitation mode is composed of several QRPA eigenmodes, 
one is generated by a u(2s~* 2 2p 3/ / 2 ) transition dominantly, and the other mostly by a i/(2s^ / 1 2 2p 1 / 2 ) 
transition. The neutron excitations take place outside of the nuclear surface reflecting the spa- 
tially extended structure of the 2s\/2 wave function. In 30 Ne, the deformation splitting of the giant 
resonance is large, and the low-lying resonance is overlapping with the giant resonance. 

PACS numbers: 21.10.Re: 21.60.Ev: 21.60.Jz 



I. INTRODUCTION 

The study of nuclei far off stability is one of the most 
active research fields in nuclear physics HSU, an< ^ 
exploring the collective motions unique in unstable nu- 
clei is one of the main issues experimentally and theo- 
retically In neutron-rich nuclei, because of the ab- 
sence of the Coulomb barrier the surface structure is 
quite different from stable nuclei. One of the unique 
structures is the neutron skin [f| [f|. Since the collec- 
tive excitations are sensitive to the surface structure, one 
can expect new kinds of exotic excitation modes asso- 
ciated with the neutron skin to appear in neutron-rich 
nuclei. One of the examples is the soft dipole excita- 
tion tii which is observed not only in light halo nu- 
clei but also in heavier 
systems jla. tl9l l20| ■ where an appreciable El strength is 
observed above the neutron threshold exhausting several 
percents of the energy- weighted sum rule (EWSR). 

The structure of the low-lying dipole state and its 
collectivity has been studied in the framework of the 
mean-field calculations by man y g roups pit . [22I [23l . [24], 

in m E3, m, m, m mi, m 113. a wiying dipole 

state in neutron-rich 26 Ne was first predicted by using 
the relativistic quasiparticle-random-phase approxima- 
tion (QRPA) in Ref. [Hj], and recently it was observed 
at RIKEN around 9 MeV, exhausting about 5% of the 
Thomas- Reiche-Kuhn (TRK) dipole sum rule [34| . In 
Ref. [13|, the QRPA was solved in the response function 
formalism. This method can treat the excitations to the 
continuum exactly by employing the Green's functions 
satisfying the out-going-wave boundary conditions, but 
an additional procedure is required to obtain the micro- 
scopic structure of the excitation mode [3o| | . 

In the present paper, we investigate the microscopic 
structure of the low-lying dipole resonance in neutron- 
rich Ne isotopes, and we discuss the isotopic dependence 
with special attention to the deformation effects. To 
this end, we have developed a deformed QRPA code in 



the matrix formulation based on the coordinate-space 
Skyrme-Hartree-Fock-Bogoliubov (HFB) theory. 

The paper is organized as follows: In Sec. [H] we ex- 
plain our method. In Sec. IIII1 we check the results of 
our new calculation scheme by comparing the existing 
QRPA results. In Sec. IIV) we present the results of the 
deformed QRPA and we discuss the microscopic struc- 
ture of the low-lying dipole state in 26 > 28 > 30 Ne. Finally, 
we summarize the paper in Sec. [Vj 



II. MODEL 

We briefly summarize here our approach (see Ref. [36| 
for details). In order to discuss simultaneously effects of 
nuclear deformation and pairing correlations including 
the continuum, we solve the HFB equations [3?], HH 
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directly in the cylindrical coordinates assuming axial and 
reflection symmetries. Here, r = v (neutron) and 7r (pro- 
ton), and r = (p, z,</>). For the mean-field Hamiltonian 
h, we employ the SkM* interaction [3^|. Details for ex- 
pressing the densities and currents in the cylindrical coor- 
dinate representation can be found in Refs. [13, SlJ. The 
pairing field is treated by using the density-dependent 
contact interaction [42, W3| , 



v pp (r,r') = V 
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with V = -390 MeV -fm 2 and g = 0.16 fin" 3 , 7=1. 
Here, g is (r) denotes the isoscalar density and P a the 
spin exchange operator. The pairing strength Vq is deter- 
mined so as to approximately reproduce the experimental 
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pairing gap of 1.25 MeV in 28 Ne obtained by the three- 
point formula [44| . Because the time-reversal symmetry 
and reflection symmetry with respect to the x — y plane 
are assumed, we have only to solve for positive f2 and 
positive z. We use the lattice mesh size Ap = Az = 0.6 
fm and the box boundary condition at p max = 9.9 fm and 
Zmax = 9.6 fm. The quasiparticle energy is cut off at 60 
MeV and the quasiparticle states up to f2' r = 1^/2^ are 
included. 

Using the quasiparticle basis obtained by solving the 
HFB equation |T]), we solve the QRPA equation in the 
matrix formulation [45j 
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The residual interaction in the particle-particle (p-p) 
channel appearing in the QRPA matrices A and B is the 
density-dependent contact interaction ((2|). On the other 
hand, for the residual interaction in the particle-hole (p- 
h) channel, we employ the Landau-Migdal (LM) approx- 
imatio n 14611 applied to the density-dependent Skyrme 
forces [471 US , 
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Here, No is the density of states and the Landau parame- 
ters are deduced from the same Skyrme force which gen- 
erates the mean field. Because the full self-consistency 
between the static mean-field calculation and the dynam- 
ical QRPA calculation is broken, we have to renormal- 
ize the residual interaction in the particle-hole channel 
by an overall factor f p h to get the spurious = 0~ 
or 1" modes (representing the center-of-mass motion) 
at zero energy (v p h — > f p h ■ v p h)- We cut the two- 
quasiparticle space at E a + Ep < 60 MeV due to the 
excessively demanding computer memory as well as the 
calculation time if we used a model space consistent with 
that adopted in the HFB calculation. Accordingly, we 
need another factor f pp for the particle-particle channel. 
We determine this factor such that the spurious = + 
mode associated with the particle number fluctuation ap- 
pears at zero energy {v pp — > /, 
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III. CHECK OF THE CALCULATION SCHEME 

In this section, we compare our results with those of 
Ref. . In this reference, the SLy4 interaction [5(| 
for the mean field and the surface-type delta interaction 
with 7 = 1.5 and Vb = —415.73 MeV-fm 3 for the pairing 
field were employed for the HFB calculation, and the 
quasiparticle energy was cut off at 50 MeV. Therefore, 
we adopt these parameters for the comparisons in this 
section. The differences between the present calculation 
and that in Ref. (4|| are the mesh size, the boundary 
condition, the cutoff energy for the QRPA calculation, 
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FIG. 1: Response function for the isoscalar quadrupole op- 
erator in 22 O. The transition strengths are smeared by us- 
ing a Lorentzian function with a width of T — 0.5 MeV. 
The renormalization factors for the QRPA calculation are 
f ph = 0.982, f pp = 1.18. The cutoff energy is 60 MeV. 



and the treatment of the spin-dependent interaction (Go 
and G ) in Eq. (fj|. In the present calculation, the spin 
transition density is treated exactly. 

In Fig. [2 we show the isoscalar quadrupole response 
function in 22 0. The first 2+ state is located at 2.8 MeV 
with B(E2 T) = 18.9 e 2 fm 4 . The experimental val- 
ues are £(2+Wp = 3.2 MeV and B(E2) c ^ p = 21 ± 8 
e 2 fm 4 [5lL I52l.l53|. In Ref. [49[, the energy and the tran- 
sition strength are E{2^) = 1.9 MeV and B{E2) = 22 
e 2 fm 4 . The energy and the transition strength of the 
low-lying collective state is quite sensitive to the cutoff 
energy for the RPA calculation 54] . In Fig. O the cut- 
off energy dependence of the renormalization factors and 
the B{E2 |) value for the 2^ state in 22 are shown. 
Even with the cutoff energy of 70 MeV, the transition 
strength for the low- lying state does not converge yet. In 
this case, the dimension of the QRPA matrix in Eq. ([3]) 
is 11726 for the = + channel and the memory size 
is 13 GB, and the CPU time is about 70,000s per each 
iteration for determining the renormalization factor f pp . 
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FIG. 2: Cutoff energy dependence of the renormalization fac- 
tors and the B(E2 f) value for the first 2 + state in 22 0. 
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If we could perform the QRPA calculation including all 
quasiparticle states obtained in the HFB calculation, the 
renormalization factor for the pairing channel f pp would 
be 1, because the p-p channel is treated self-consistently 
between the HFB and the QRPA calculations. 

The peak position of the giant resonance is located 
slightly higher than in Ref.[49(. The non-collective two- 
quasiparticle states around 6 and 7 MeV are consistent 
between the two calculations. The energy- weighted sum 
(1.867 x 10 4 MeVfm 4 ) overestimates by about 13.9% the 
EWSR value (1.638 x 10 4 MeVfm 4 ). The overshooting of 
the EWSR for the isoscalar quadrupole mode in the LM 
approximation was pointed out in Ref [55| . 



dipole operator used in the present calculation is 



Fik = e 
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and the response functions are calculated as 

s(£)=EE r/2 l(z|/w|0>|2 
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IV. RESULTS AND DISCUSSION 

We now discuss the properties of 26 > 28 > 30 Ne nuclei cal- 
culated with the SkM* interaction. We summarize in 
Table U the ground state properties of these Ne iso- 
topes obtained by solving Eq. (Q}. The ground state 
is slightly deformed in 26 Ne and 28 Ne, and we obtain 
a well-deformed ground state for 30 Ne. The values in 
parentheses are experimental pairing gaps extracted by 
the three-point mass difference formula [4J| using the ex- 
perimental binding energies taken from Ref. [561 ]. We 
define the deformation parameter @2 and average pairing 

gap (a) (Hmmiii as 



We can clearly see a resonance structure at around the 
excitation energy of 8-9 MeV, together with the giant res- 
onance at 15 — 20 MeV. Because of the small deformation 
the K splitting is small and smeared out. 

In Fig. 31 we show the transition strengths in the low- 
energy region. The neutron emission threshold is 6.35 
MeV, and the resonance which is composed of several 
discrete states appears just above the threshold. In con- 
trast to the low- lying quadrupole state in 22 O, the tran- 
sition strengths for the dipole states in this region con- 
verge at the cutoff energy of about 40 MeV. We made a 
detailed analysis of the QRPA eigenmodes and show in 
Tables IIIIIIII the microscopic structures of the IC = 0~ 
state at 8.25 MeV and the = 1" state at 8.76 MeV, 
which have the largest transition strength for each sec- 
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where g(r) is the pairing density. 

Fig. [3] shows the response functions for the isovector 
dipole mode in neutron-rich Ne isotopes. The isovector 



TABLE I: Ground state properties of 



Ne obtained by 



the deformed HFB calculation with the SkM* interaction and 
the surface-type pairing interaction. Chemical potentials, de- 
formations, average pairing gaps and root-mean-square radii 
for neutrons and protons are listed. 





2b Ne 


2S Ne 


3U Ne 


A„ (MeV) 


-4.60 


-3.06 


-2.90 


K (MeV) 


-14.8 


-17.0 


-19.9 


ft 


0.08 


0.12 


0.32 




0.14 


0.20 


0.39 


(A„) (MeV) 0.0 (0.70) 1.27 (1.24) 1.34 (1.30) 


<A^) (MeV) 


1.04 


0.87 


0.0 


V(r 2 )» (fm) 


3.20 


3.35 


3.53 


V(r 2 )« (M 


2.93 


2.98 


3.08 



TABLE H: QRPA amplitudes for the A" 1 = 0~ state in 
26 Ne at 8.25 MeV. This mode has the proton strength 

-2 „2r m 2 



e fm , the neutron strength B{Q V 1) 
e 2 fm 2 , and the isovector strength B(Q IV 1 



x 10 

,2f m 2 



B{E1) = 2.£ 
2.89 x 10~ 2 

1.17 x 10-Vfm 2 , and the sum of backward-going amplitude 
I]|V a/3 | 2 = 4.33 x 10" 3 . The single- (quasi) particle levels 
are labeled with the asymptotic quantum numbers [Nn^K\Q.. 
Only components with A 2 ^ — Y^g > 0.001 are listed. Two- 
quasiparticle excitation energies are given by E a + E/s in MeV 
and two-quasiparticle transition matrix elements Qio, a /3 in 
e-fm. In the row (i), the label z/1/2 - denotes a non-resonant 
discretized continuum state of neutron Q n = l/2~ level. 









E a + Eg 




Ql0,af3 




a 


P 


(MeV) 


y-2 V 2 


(e- fm) 


(a) 


v[3W] 1/2 


i/[211]l/2 


8.15 


0.670 


-0.309 


(b) 


v [330] 1/2 


i/[220]l/2 


11.4 


0.020 


-0.397 


(c) 


j/[312]5/2 


v [202] 5/2 


11.2 


0.006 


-0.239 


(d) 


i/[321]3/2 


i/[211]3/2 


11.3 


0.006 


0.338 


(e) 


i/[330]l/2 


i/[211]l/2 


6.54 


0.003 


-0.118 


(f) 


i/[312]3/2 


i/[211]3/2 


12.8 


0.002 


-0.014 


(g) 


i/[301]l/2 


i/[211]l/2 


9.32 


0.002 


-0.117 


(h) 


i/[200]l/2 


i/[101]l/2 


14.0 


0.002 


-0.241 


(i) 


1/1/2" 


i/[211]l/2 


12.6 


0.002 


-0.068 


(j) 


n [220] 1/2 


tt[101]1/2 


7.96 


0.265 


0.0085 


(k) 


tt[330]1/2 


tt [220] 1/2 


13.4 


0.008 


-0.329 


(1) 


tv [220] 1/2 


tt[110]1/2 


14.1 


0.008 


-0.346 
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FIG. 3: Response functions for the isovector dipole operator in 26 ' 28 ' 30 N e . The dotted, dashed and solid lines correspond to the 
K 71 = CP, K 7 ' = 1~ and total responses, respectively. For the K 7 ' = 1~ response, the transition strengths for the K 7 ' — ±1 _ 
states are summed up. The transition strengths are smeared by using F = 1 MeV. The renormalization factors for the QRPA 
calculation are f ph = 0.919, 0.880 and 0.929 for 26 . 28 ' 30 Ne and f pp = 1.225 for all nuclei. 
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FIG. 4: Isovector dipole transition strengths in 26 Ne for the 
K 71 — 0~ (the upper) and JC = 1 _ (the lower) states. Un- 
derlying discrete states are shown together with the smeared 
response functions. The arrow indicates the neutron emission 
threshold E th = 6.58 MeV. 



tor. In the Tables, single-(quasi)particle states are la- 
beled with the asymptotic quantum numbers [TV 713 A] SI 
just for convenience. It should be noted that the asymp- 
totic quantum numbers are not good quantum numbers 
because the deformation is not so large. 

For the K* = 0~ state at 8.25 MeV, the dominant 
component is the i/[211]l/2 — » i/[310] 1/2 transition, cor- 
responding to i / (2s7, 2 2p3/ 2 ). Particle-hole excitations of 

(b), (c) and (d) correspond to v{ld§ ^1/7/2) excitation, 
which have 3.2% contribution in total. The rows (e), (f), 
(g) and (h) correspond to v(2s~J 2 l} 7 /2 ), v{ld~j 2 2p 3/2 ), 

!/(2s7 1 2pi/2) and v{\p~j 2 ld 3 / 2 ) excitations, respectively. 
Two-quasiparticle proton excitations, furthermore, have 
an appreciable contribution; the rows (j), (k), and (1) 
correspond to the hole-hole like ir{lpi/2 <8> 1^5/2) excita- 
tion, particle-hole like Tr(ld~7} 2 lf 7 / 2 ), and Tr(ld~} 2 lp 3 / 2 ) 



TABLE III: Same as Table [TT] but for the K 7 ' = 1" state in 
26 Ne at 8.76 MeV. This mode has B(E1) = 1.65 x 10" 2 e 2 fm 2 , 
B{Q"l) = 3.58 x 10 _2 e 2 fm 2 , B(Q IV 1) = 1.00 x 10 _1 e 2 fm 2 , 
and J2 \ Y ap\ 2 = 2 - 93 x 10~ 3 . 









E a + Ep 




Qll,af3 




a 




(MeV) 


V 2 V 2 


(e- fm) 


(a) 


z/[312]3/2 


i/[211]1/2 


8.68 


0.849 


0.339 


(b) 


i/[310]l/2 


v [211] 1/2 


8.16 


0.040 


-0.131 


(c) 


i/[301]l/2 


v [211] 1/2 


9.32 


0.010 


0.294 


(d) 


i/[321]3/2 


i/[220]1/2 


12.0 


0.007 


0.250 


(e) 


^[303] 7/2 


^[202] 5/2 


12.1 


0.006 


0.414 


(f) 


v [330] 1/2 


v [220] 1/2 


11.4 


0.004 


-0.127 


(g) 


^[312]5/2 


i/[211]3/2 


12.1 


0.004 


0.348 


00 


i/[321]3/2 


^[202]5/2 


10.3 


0.001 


-0.010 


(i) 


i/[321]3/2 


i/[202]5/2 


11.8 


0.003 


-0.214 


(J) 


1/ [330] 1/2 


i/[211]1/2 


6.54 


0.003 


-0.081 


(k) 


f[321]3/2 


i/[211]1/2 


7.14 


0.001 


0.106 


(1) 


7T [220] 1/2 


7r[101]l/2 


7.96 


0.037 


0.0095 


(m) 


tt[211]3/2 


tt[101]1/2 


7.95 


0.015 


-0.011 


(n) 


tt[321]3/2 


tt[220]1/2 


14.0 


0.004 


0.313 


(0) 


tt[312]5/2 


tt[211]3/2 


14.7 


0.002 


-0.338 


(p) 


tt[211]3/2 


7r[110]l/2 


14.1 


0.002 


0.280 


(q) 


7t[211]1/2 


tt[101]1/2 


11.5 


0.002 


-0.256 



excitations, respectively. The K 71 = 1~ state has a sim- 
ilar structure to the K n = 0~ state. The main com- 
ponent is v(2s~j~ 2 2p 3 / 2 ), which corresponds to (a) and 
(b) in Table HTT1 With a small contribution, many other 
neutron particle-hole and proton two-quasiparticle exci- 
tations build the excitation mode at 8.76 MeV. The reso- 
nance is also composed of the K 71 = 1~ mode appearing 
at 9.40 MeV. This mode is dominantly (97.6%) generated 
by the ^[211]l/2 — > z/[301]l/2 transition corresponding to 
the v(2s~^J 2 2p 1 / 2 ) transition. 

Preliminary calculations of deformed QRPA using the 
Gogny interaction [6l[, and the relativistic deformed 
QRPA [62] show that the low-lying dipole state is domi- 
nantly constructed by the v{2s~^ 2 2p 3 j 2 ) configuration. 
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FIG. 5: Transition densities in 26 Ne for the K 7 ' = 0~ state at 
8.25 MeV (upper panels), and for the K n = 1~ state at 8.76 
MeV (lower panels). Solid and dotted lines indicate positive 
and negative transition densities, and the contour lines are 
plotted at intervals of 3 x 10~ 4 fm -3 . The thick solid lines 
indicate the neutron and proton half density, 0.058 fm -3 and 



0.036 fm , respectively. 



v[211]1/2->[310]1/2- 




FIG. 6: Same as Fig. [5] for the unperturbed transition density 
of i/[211]l/2 -» [310] 1/2 excitation. 



In Fig. we show transition densities for the K* = 
state at 8.25 MeV, and for the K* = 1" state at 8.76 
MeV. These transition densities are quite different from 
the classical picture of the isovector giant resonances. 
They have an isoscalar character in the surface region of 
the nucleus. On the other hand, outside of the nucleus, 
neutrons have an oscillation and the neutron excitation 
is dominant. Furthermore, the neutron excitations take 
place in the low-density region around 6 fm, namely, this 
mode has a unique picture of vibration of the neutron 
skin. The spatially extended structure of the z/[211]l/2 
state is responsible for this tail of the neutron transi- 
tion density. In order to clearly see the spatial struc- 
ture of v [2 11] 1/2, the unperturbed transition density of 
z/[211]l/2 -> [310] 1/2 is shown in Fig.© It is clear that 
the wave function extends far outside of the nucleus, and 
the extension is larger than in the dipole state of Fig. [5] 
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FIG. 7: Energy weighted sum of the isovector dipole strength 
function. The solid and the dotted lines are calculations with 
the energy cutoff at 60 and 30 MeV. The horizontal lines 
show the classical TRK, and the RPA sum rule including the 
enhancement factor, mi = mi (1 + ft). Here k = 0.32 for the 
SkM* interaction in 26 Ne. 



Furthermore, the structure around the surface is differ- 
ent between the unperturbed one and that obtained in 
the QRPA. These differences are generated by the QRPA 
correlations; the low-lying dipole state possesses a collec- 
tive nature, which is small but finite. 

Fig. shows the energy weighted sum of the isovec- 
tor dipole strength function together with the sum rule 
values represented by the horizontal lines. The calcu- 
lated sum satisfies 89.2% of the EWSR value including 



the enhancement factor 



k: mi 



'(1 + k) [26J. The 



enhancement factor comes from the momentum depen- 
dence of the Skyrme density functionals. The effect of 
the explicit treatment of the momentum dependence for 
the EWSR was discussed in the discretized-continuum 



QRPA [63J and the continuum QRPA 55] for the spheri- 
cal systems. In the present calculation, we treat the mo- 
mentum dependence in the LM approximation. There- 
fore, discrepancy between the calculation and the EWSR 
value comes from this treatment of the momentum de- 
pendence. This point remains to be improved, and it is 
discussed in Ref. [rJi ]. 

In the present calculation, the energy-weighted sum up 
to 10 MeV is 5.51 MeV-e 2 fm 2 , corresponding to 6.0% of 
the TRK sum-rule value, 4.6% of the EWSR including 
the enhancement factor and 5.1% of the calculated sum. 
These values are consistent with the experiment [34l ]. In 
Fig. we also show the energy-weighted sum calculated 
with the energy cutoff at 30 MeV (dotted line). In the 
giant resonance region, two calculations give different re- 
sults, while they are almost identical in the low-energy re- 
gion. This is because the collectivity of the low-lying res- 
onance is small, and consequently the transition strength 
is not very sensitive to the cutoff energy. 

Before going to the neighboring nuclei, it should 
be noted that we obtain the collective octupole 
state at about 5.2 MeV, below the neutron thresh- 
old, with B(E3 t) = 2458e 2 fm 6 , which corresponds 
to about 61 in Weisskopf units and the isoscalar 
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FIG. 8: Same as Fig. g] but for 28 Ne and 30 Ne. The arrows 
indicate the neutron emission threshold E t h = 4.44 and 4.51 
MeV. 



transition strength is 2.60 x 10 4 frn 6 . The lowest 
K* = 0~ state is located at 5.03 MeV and the sum 
of the backward-going amplitudes is 0.099. This 
state is generated by i/[211]l/2 -> [330] 1/2(53.2%), 
z/[202]5/2 -> [312]5/2(6.8%), i/[211]l/2 -> 
[310]l/2(6.5%), i/ [220] 1/2 -> [330]l/2(3.7%), and 
7r[101]2/l®[220]l/2(14.8%),7r[101]3/2(8)[211]3/2(2.2%), 
tt[220]1/2 (g) [330]1/2(2.1%), tt[101]1/2 <g> [211]1/2(1.0%). 



TABLE IV: Same as Table [TT] but for the IC = 0" state at 
8.14 MeV in 28 Ne. This mode has B(E1) = 2.62 x 1(T 2 e 2 fm 2 , 
B(Q"1) = 2.29 x 10" 2 e 2 fm 2 , B(Q IV 1) = 9.80 x 10~ 2 e 2 fm 2 , 
and J2 \ Y c*0\ 2 = 9.77 x 10~ 3 . In the rows (k), (1), (m), (n) and 
(t), the labels vl/2~ , vl/2 + and 7iT/2 + denote non-resonant 
discretized continuum states of neutron tt n = 1/2 - and l/2 + 
levels and proton l/2 + level. 
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-0.463 
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-0.441 
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B. 28Ne and 30Ne 

The central panel in Fig.[3]shows the response function 
in 28 Ne. In the low-energy region, we can see a two-bump 
structure at around 7 and 8 MeV. Because the deforma- 
tion is small as in 26 Ne, we cannot see a splitting of the 
giant resonance. In Fig. [5] we show the low-energy part 
of the strength functions. In the = 0~ states, there 
is a prominent peak at 8.1 MeV with a strength of 0.098 
e 2 fm 2 . The strength distribution is fragmented for the 
= 1" mode, but correspondingly, we can see an eigen- 
mode at 8.9 MeV with the largest transition strength of 
0.058 e 2 fm 2 . 

We show in Table |IV] the QRPA amplitude for the 
K n = 0~ state at 8.14 MeV in 28 Ne. The main com- 
ponent is the neutron two-quasiparticle excitation of 
i/([310]l/2 (g> [211]l/2) corresponding to v(2s~} 2 2p 3/2 ). 

Two quasiparticle excitations of (b) and (c) in Table ITVl 
correspond to u(ld^} 2 lf 7 / 2 ), and (d): v(2s~j 2 lf 7 / 2 ), (e) 

and (f): v(ld~j 2 2p 3/2 ), (g): v{\d~j 2 \f 7 / 2 ), (h.) and (i): 
i/(l(ir, 2 2p3/2), and (j): u(ld^, 2 2pi/ 2 ) excitations, respec- 
tively. The proton excitation of Tt(lpi/ 2 (S> W5/2) has an 
appreciable contribution as in 26 Ne. 

The lower energy resonance at around 7 MeV is de- 
scribed by three eigenstates as shown in Fig. [5J The 
lower state at 6.70 MeV, which has an isovector strength 



of 0.028 e 2 fm 2 , is mainly generated by v(ld 3 f 2 2p 3 / 2 ) 
(87.5%), 1/(25-^1/7/2) (5.7%) and u(ld-^2p 1/2 ) (2.4%). 
The state at 6.96 MeV with B(Q w l) = 0.044 e 2 fm 2 is 
almost a single p-h excitation of u(ld^} 2 2pi/ 2 ) (90.2%), 
and the state at 7.19 MeV with 0.021 e 2 fm 2 is gen- 
erated dominantly by proton h-h like excitation of 
7r(lpi/ 2 <8 ld 5/2 ) (57.8%), together with v(2s^ 2 lf 7 / 2 ) 

(18.6%), 1/(1^2^3/3) (6.4%), ^(lrf-/ 1 ^/^) (5.2%), 

v(2s-f 2 2p 3/2 ) (4.5%) and v(ld-f 2 2p 1/2 ) (1.8%). 

Therefore, the higher-energy resonance at 8 MeV has 
a similar structure to that in 26 Ne; v{2s^, 2 2p 3 j 2 ) and 

7r (lpi/2 ® 1^5/2) excitations are dominant, and the lower- 
energy resonance is generated by different eigenmodes. 

The K* = I" state at 8.9 MeV in 28 Ne has a similar 
structure to that in 26 Ne and K v = 0" state in 28 Ne, cor- 
responding mainly to the neutron two-quasiparticle exci- 
tation of u(2sy 2 2p 3 / 2 ), with 64.0% contribution. In ad- 
dition to this neutron p-h like excitation, the following ex- 
citations have an appreciable contribution; u(ld^, 2 lf 7 / 2 ) 

(13.9%), v{2s-f 2 lf 7/2 ) (2.8%), 7r(lp^ld 5/2 ) (7.9%) and 
7r(lp 3/2 <8> ld 5/2 ) (1.4%). 

Finally, we discuss the dipole state in 30 Ne. Compared 
to the response functions in 26 Ne and 28 Ne, that for 30 Ne 
is quite different, because this nucleus is well deformed 
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TABLE V: Same as Table HQ but for the IC = 1" state at 
6.69 MeV in 30 Ne, and only components with X^g - Y£ p > 
0.01 are listed. This mode has B(E1) = 1.58 x 10~ 2 e 2 fm 2 , 



B{Q v l) = 3.98 x 10 _; Vn 
and \Y a g\ 2 = 8.82 x 10" 



B(Q 1V 1) = 1.06 x 10~Vfm 2 
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FIG. 9: Same as Fig.[5]but for the K* = 0" state at 8.1 MeV 
and for the K n = 1" state at 6.7 MeV in 30 Ne. 



as shown in Table HI The giant resonance is split into 
K 77 = CP and 1~ mode, and the split giant resonance has 
an overlap with the low-lying resonance below 10 MeV. In 
the right panel of Fig. [5J we show the strength distribu- 
tion below 10 MeV in 30 Ne. For the K"* = 0~ mode, we 
can see a prominent peak at 8.1 MeV possessing a large 
isovector El strength of 0.48 e 2 fm 2 . This state is mainly 
generated by i/[211]l/2 -» [310]l/2 (38.4%) and the neu- 
tron excitation from [330] 1/2 to the non- resonance con- 
tinuum state (37.9%), together with the proton excita- 
tion of tt[330]1/2 -> [220] 1/2 (6.1%). 

In Fig. [51 the transition density for the K v = 0~ state 
is shown. The transition density of protons are quite sim- 
ilar to that in Fig. [5] For the neutrons, we can easily see 
the effect of mixing of the excitation into the continuum 
state; the transition density has large spatial extension. 
Furthermore, comparing to Fig. O this K* = 0" state 
still possesses a structure similar to the low-lying dipole 
state in 26 Ne. 

For the K 71 — l~ state, we can see a prominent peak at 
6.69 MeV possessing an isovector strength of 0.11 e 2 fm 2 . 
This state has a different structure to the dipole states 
discussed above. In Table El we show its microscopic 
structure. This state has a collective nature in a sense 
that a number of two-quasiparticle excitations have an 
appreciable contribution; in the present case, eight of the 
neutron excitations have a contribution larger than 1%. 
In the lower panel of Fig. [51 the transition density of 
this state is shown. This mode has also a characteristic 
feature that the neutron and proton contribution have 
an isoscalar nature around the surface region, and the 
neutron excitation is dominant outside of the nucleus. 

It is difficult to link directly with the low-lying dipole 
states in 26 Ne or 28 Ne, because the deformations are quite 



different in 30 Ne and in the other two nuclei. The main 
component of this K v — 1~ state is (a):v[200]l/2 — > 
[312J3/2 and (b):i/[200]l/2 -» [310]l/2. These p-h exci- 
tations are ^(ld^^p^^) m the spherical limit. In this 
sense, the lower-energy resonance in 28 Ne is connected to 
this collective K 77 = 1~ state in 30 Ne. 

V. SUMMARY 

We have investigated a new framework of the deformed 
QRPA based on the Skyrme density functionals and the 
Landau-Migdal approximation. With this method, we 
have made a detailed analysis of the low-lying dipole 
states in neutron-rich 26 > 28 > 30 Ne. In these nuclei, we ob- 
tain the excitation mode at 8 — 8.5 MeV. The low-lying 
resonance is composed of several QRPA eigenmodes. In 
26 Ne, not only the v(2s^ 2 2p 3 / 2 ) transition but also the 

v{2s~J 2 2pi/2) transition contribute to generating the res- 
onance. In 28 Ne and 30 Ne, the i/[211]l/2 -> [310] 1/2 
excitation still plays a major role. Each eigenmode is, 
however, not purely a single particle-hole excitation, it 
has a small contribution of the other neutron excitations 
and proton excitations as well. 

We have clearly shown the spatially extended structure 
of the f[211]l/2 (2si/ 2 ) state, and that it is responsible 
for the oscillation of transition density of neutrons out- 
side of the nucleus. In the well deformed nucleus 30 Ne, 
the deformation splitting of the giant resonance is large 
and the low-lying resonance overlaps with the giant res- 
onance. For the — 1~ state, we furthermore obtain 
a collective dipole mode at 6.7 MeV. 
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